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The confusion and ambiguity encountered by students, in understanding virtual displacement 
and virtual work, is addressed in this article. A definition of virtual displacement is presented that 
allows one to express them explicitly for both time independent and time dependent constraints. It 
is observed that for time independent constraints the virtual displacements are the displacements 
allowed by the constraints. However this is not so for a general time dependent case. For simple 
physical systems, it is shown that, the work done on virtual displacements by the constraint forces 
is zero in both the situations. For allowed displacements however, this is not always true. It is also 
demonstrated that when constraint forces do zero work on virtual displacement, as defined here, we 
have a solvable mechanical problem. We identify this special class of constraints, physically realized 
and solvable, as the ideal constraints. The concept of virtual displacement and the principle of zero 
virtual work by constraint forces are central to both Lagrange's method of undetermined multipliers, 
and Lagrange's equations in generalized coordinates. 



I. INTRODUCTION 

Almost all graduate level courses in classi- 
cal mechanics include a discussion of virtual 
displacemen1ii*2i2i!i£i£iL2i2ii2iii and Lagrangian 
dynamicob 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ^ 1 ! 12 ^ 3 . From the concept 
of zero work by virtual displacement the Lagrange's 
equations of motion are derived. However, the definition 
of virtual displacement is rarely made precise and often 
seems vague and ambiguous to students. In this article 
we attempt a more systematic and precise definition, 
which not only gives one a qualitative idea of virtual 
displacement, but also allows one to quantitatively 
express the same for any given constrained system. 
We show that in a number of natural systems, e.g., 
particle moving on a frictionless slope, pendulum with 
moving point of suspension, the work done by the forces 
of constraint due to virtual displacement is zero. We 
also demostrate that this condition is necessary for the 
solvability of a mechanical problem. Hence we propose 
such systems as an important class of natural systems. 



A. Ambiguity in virtual displacement 

In the following we try to classify the difficulties faced 
by a student in understanding the definition of virtual 
displacement. 

1. It is claimed that (i)a virtual displacement 5r is 
consistent with the forces and constraints imposed 
on the system at a given instant tX>; (ii) a virtual dis- 
placement is an arbitrary, instantaneous, infinites- 
imal change of position of the system compatible 
with the conditions of constraint; (iii) virtual dis- 
placements are, by definition, arbitrary displace- 



ments of the components of the system, satisfy- 
ing the constraint 3 ; (iv) virtual displacement does 
not violate the constraints®; (v) we define a vir- 
tual displacement as one which does not violate the 
kinematic relations^; (vi) the virtual displacements 
obey the constraint on the motion 11 . These state- 
ments imply that the virtual displacements satisfy 
the constraint conditions, i.e., the constraint equa- 
tions. However this is true only for time indepen- 
dent (sclerenomous) constraints. We shall show 
that for time dependent (rheonomous) constraints, 
such as a pendulum with moving support, this def- 
inition would violate the zero virtual work condi- 
tion. 

2. It is also stated that (i) virtual displacement is to 
be distinguished from an actual displacement of the 
system occurring in a time interval dim; (ii) it is 
an arbitrary, instantaneous, change of position of 
the system^; (iii) virtual displacement Sr takes place 
without any passage of time®, (iv) virtual displace- 
ment has no connection with the time - in contrast 
to a displacement which occurs during actual mo- 
tion, and which represents a portion of the actual 
path 9 ; (v) one of the requirements on acceptable 
virtual displacement is that the time is held fixed—. 
We even notice equation like l Sxi — dxi for dt — 0' 8 . 
The above statements are puzzling to a student. If 
position is a continuous function of time, a change 
in position during zero time has to be zero. In 
other words, this definition implies that the vir- 
tual displacement cannot possibly be an infinitesi- 
mal (or differential) of any continuous function of 
time. In words of Arthur Haas: since its (virtual 
displacement) components are thus not fucntions of 
the time, we are not able to regard them as differ- 
entials, as we do for the components of the element 
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of the actual path—. We shall show that virtual dis- 
placement can be looked upon as a differential, it 
is indeed a differential increment in virtual velocity 
over a time dt, Eq.©. 

3. It is also stated that (i) virtual displacements do 
not necessarily conform to the constraints-; (ii) the 
virtual displacements 5q have nothing to do with 
actual motion. They are introduced, so to speak, as 
test quantities, whose function it is to make the sys- 
tem reveal something about its internal connections 
and about the forces acting on jja; (iii) the word 
"virtual" is used to signify that the displacements 
are arbitrary, in the sense that they need not corre- 
spond to any actual motion executed by the systerA; 
(iv) it is not necessary that it (virtual displacement) 
represents any actual motion of the system 6 ; (v) 
it is not intended to say that such a displacement 
(virtual) occurs during the motion of the particle 
considered, or even that it could occw&; (vi) virtual 
displacement is any arbitrary infinitesimal displace- 
ment not necessarily along the constrained path 5 . 
From the above we understand that the virtual 
displacemnts do not always satisfy the constraint 
equations, and they need not be the ones actually 
realized. We shall see that these statements are 
consistent with physical situations, but they cannot 
serve as a satisfactory definition of virtual displace- 
ment. Statements like: "not necessarily conform to 
the constraints" or "not necessarily along the con- 
strained path" only tell us what virtual displace- 
ment is not, they do not tell us what it really is. 
Reader should note that there is a conflict between 
the claims under items 1 and 3. It is not clear from 
the above, whether the virtual displacements sat- 
isfy the constraints, i.e., the constraint equations 
or not. 

4. Virtual displacement is variously described as: ar- 
bitrary, virtual, and imaginary 1 - 2 - 3 5 6 . These adjec- 
tives make the definition somewhat mysterious to 
a student. 

Together with the above ambiguities, students are of- 
ten confused as to whether it suffices to understand vir- 
tual displacement as an abstract concept, or they need to 
have a quantitative definition. Some students appreciate 
that the virtual displacement as a vector should not be 
ambiguous. The principle of zero virtual work is used 
to derive Lagrange's equations. For a particle under con- 
straint this means that the virtual displacement is always 
orthogonal to the force of constraint. 

At this stage a student may get further puzzled. 
Should he take the forces of constraint as supplied, and 
then the principle of zero virtual work as a definition 
of virtual displacement ? In that case the principle re- 
duces merely to a definition of a new concept, namely 
virtual displacement. Or should the virtual displacement 
be defined from the constraint conditions independently 



? The principle of zero virtual work may then be used to 
calculate the forces of constraint that ensure constraint 
condition throughout the motion. 

II. VIRTUAL DISPLACEMENT AND FORCES 
OF CONSTRAINT 

A. Constraints and Virtual displacement 

Let us consider a system of constraints that are ex- 
pressible as equations involving positions and time. They 
represent some geometric restrictions (holonomic) either 
independent of time (sclerenomous) or explicitly depen- 
dent on it (rheonomous). Hence for a system of N par- 
ticles moving in three dimensions, a system of (s) holo- 
nomic, rheonomous constraints are represented by func- 
tions of Yk and (t), 

f l {v 1 ,Y 2 ,...,Y N ,t) = Q, i = 1, 2, ...,s (1) 

Each constraint of this form imposes a restriction on the 
possible or allowed velocities, which must satisfy, 

=»• • (2 » 

It is worth noting at this stage that there are many, in 
fact infinitely many, allowed velocities, since we have im- 
posed only (s) number of constraints on (3N) scalar com- 
ponents of the allowed velocity vectors. An infinitesimal 
displacement over time (dt) due to allowed velocities will 
be called the allowed infinitesimal displacement or simply 
allowed displacement. 

dr k =v k dt k = l,2,...,N (3) 

Allowed displacements dr k together with differential of 
time (dt) satisfy constraint equations similar to Eq.(|2")l. 

k— 1 v 7 

As there are many allowed velocities we have many al- 
lowed infinitesimal displacements. We propose to define 
virtual displacement as the difference between any two 
such (unequal) allowed displacements, 

Sr k =dr k -dr' k , k = l,2,...,N (5) 

This definition is motivated by the possibility of identi- 
fying a special class of 1 ideal constraints^ (sec. He), and 
verifying : the principle of zero virtual work 1 in common 
physical examples (sec. III). It may be noted that, by 
this definition, virtual displacement 5r k is not a change 
in position in zero time. It is rather the difference of any 
two allowed displacements during a time dt. 

5r k = (v k -v' k )dt, k = l,2,...,N (6) 
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The difference of two allowed velocities v k = — vl may 
be defined as the virtual velocity. 

The virtual displacements thus defined satisfy the ho- 
mogeneous part of the constraint equation Eq.Q) (i.e., 
with dfi/dt = 0). 

N „, 

J2^-5r k = 0, » = l,2,..., a (7) 

fc=i rfe 

The absence of the (dfi/dt) in the above equation, 
Eq.Q, gives the precise meaning to the statement that 
virtual displacements are the allowed displacements in the 
case of frozen constraints. The constraints are frozen in 
time in the sense that we make the (dfi/dt) term zero, 
though the dfi/dr k term still involves time. In the case of 
stationary constraints, i.e., / (ri, . . . , rjv) = 0, the virtual 
displacements are identical with allowed displacements as 
(dfi/dt) is zero. 

B. Existence of forces of constraints 

In the case of an unconstrained system of TV particles 
described by position vectors (rfe) and velocity vectors 
(vfe ) , the motion is governed by Newton's Law, 

m k a k = F k (r l ,v l ,t), k,l= 1,2,..., N (8) 

where m k is the mass of the fcth particle, a k is its ac- 
celeration and F k is the total external force acting on it. 
However, for a constrained system, the equations of con- 
straint, namely Eq.JTJ, impose the following restrictions 
on the allowed accelerations, 

^ = 1,2,..., S (9) 

Given rfe, v k one is no longer free to choose all the ac- 
celerations a/e independently. Therefore in general the 
accelerations afe allowed by Eq. © are incompatible with 
Newton's Law, 

m fc a fc =F fc , k = l,2,...,N 

This implies that during the motion the constraint con- 
dition cannot be maintained by the external forces alone. 
Physically some additional forces, e.g., normal reaction 
from the surface of constraint, tension in the pendulum 
string, come into play to ensure that the constraints are 
satisfied, Hence one is compelled to introduce forces of 
constraints Rfe and modify the equations of motion as, 

m fc a fe = Ffe + R fc , k = l,2,...,N (10) 

Now the problem is to determine the motion of TV parti- 
cles, namely their positions (r k (t)), velocities (v k (t)) and 
the forces of constraints (Rfe), for a given set of external 



forces Ffe, constraint equations (/j(ri,r2, . . . , tn, t) = 0, 
i = 1,2, ... ,s) and initial conditions (rfc(O), Vfe(0)). It is 
important that the initial conditions are also compatible 
with the constraints. 

There are a total of (67V) scalar unknowns, namely the 
components of r k (t) and Rfe, connected by (37V) equa- 
tions of motion, Eq. (|10fl . and (s) equations of constraints, 
Eq.JTJ. For (67V > 37V + s) we have an under-determined 
system. Hence to solve this problem we need (37V — s) 
additional scalar relations. 



C. Solvability and ideal constraints 

In simple problems with stationary constraints, e.g., 
motion on a smooth stationary surface, we observe that 
the allowed displacements are tangential to the surface. 
The virtual displacement being a difference of two such 
allowed displacements, is also a vector tangential to it. 
The force of constraint, so called 'normal reaction', is 
perpendicular to the surface. Hence the work done by 
the constraint forces on allowed as well as virtual dis- 
placement is zero, 

JV N 

• dr k =0, J2 Rk Srk = 

fc=i fe=i 

When the constraint surface is in motion, the allowed 
velocities, and hence the allowed displacements are no 
longer tangent to the surface (see sec. III). The virtual 
displacement remains tangent to the constraint surface. 
If the forces of constraint can still be assumed normal to 
the instantaneous position of the surface, we have zero 
virtual work. However note that the work by constraint 
forces on allowed displacements is not zero. 

JV N 

^Rfe-dr fc ^0, ^R fe -(5rfe = (11) 

fe=l fe=l 

In a number of physically interesting simple problems, 
such as, motion of a pendulum with fixed and moving 
support, motion of a particle along a stationary and mov- 
ing slope, we observe that the above interesting relation 
between the force of constraint and virtual displacement 
holds (see sec. III). Out of the above 37V virtual dis- 
placements, only n = 37V — s are independent. If the (s) 
dependent quantities are expressed in terms of remaining 
n = 37V — s independent objects we get 

n 

Y,Rj-Sxi = Q (12) 

3=1 

where the independent components of rfe. Rj are 

the coefficients of Sxj, and are composed of different Rfe. 
Since the above components of virtual displacements Sxj 
are independent, one can equate each of their coefficients 
to zero (Rj = 0). This brings in (37V — s) new scalar 
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conditions or equations and the system is solvable (not 
under-determined) again. 

Thus we find a special class of constraints which is 
observed in nature (see sec. Ill) and which gives us a 
solvable system. We call this special class of constraints, 
satisfying zero virtual work principle by constraint forces, 
i.e., ^2 k Rfc • <5rfc = 0, the ideal constraint. 

Our interpretation of the principle of zero virtual work, 
as a definition of an ideal class of constraints, find sup- 
port in Sommcrfcld. In his words, "a general postulate of 
mechanics: in any mechanical systems the virtual work 
of the reactions equals zero. Far be it from us to want to 
give a general proof of this postulate, rather we regard it 
practically as definition of a mechanical system ". 

III. EXAMPLES OF VIRTUAL 
DISPLACEMENTS 

A. Simple Pendulum with stationary support 

The motion of the pendulum is confined to a plane and 
the bob moves at a fixed distance from the point of sus- 
pension. The equation of constraint by Eq.JU therefore 
is, 

f(x, y, t) = x 2 + y 2 - rl = 

Whence 




FIG. 1: Allowed and virtual displacements for a pendulum 
with stationary support 

Hence the constraint equation for allowed velocities 
(compare Eq.J2J)) is, 

x ■ v x + y ■ v y = 

Hence the allowed velocity (v x , v y ) is orthogonal to the 
instantaneous position (x,y) of the bob relative to sta- 
tionary support. The same may also be verified taking a 
plane polar coordinate. 



The allowed displacements are always collinear to al- 
lowed velocities. Virtual displacement being difference 
of two allowed displacements, is also a vector collinear 
to the allowed velocities, hence tangential to the line of 
suspension. 

dr = vdt, dr' — v'dt 

5r=(v- v')dt 

We may assume that the string of the pendulum provides 
a tension (T) but no shear (ideal string). We get zero 
work by tension due to allowed and virtual displacements, 

T • dr = 0, T • 5r = 

B. Simple Pendulum with moving support 

Let us first consider the case when the support is mov- 
ing vertically with a velocity u. The motion of the pen- 
dulum is still confined to a plane. The bob moves keeping 
a fixed distance from point of suspension. The equation 
of constraint is, 

f(x,y,t)=x 2 + (y-ut) 2 -r 2 = 

where u is the velocity of the point of suspension along 
a vertical direction. 




FIG. 2: Allowed and virtual displacements for a pendulum 
with moving support 

Whence 

^-=2x, ^f = 2(y-ut), ^- = -2u(y - ut) 
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Hence the constraint equation gives, 

x ■ v x + (y — ut) ■ v y — u(y — ut) = 

or, 

x ■ v x + (y - ut) ■ (v y — u) = 

Hence the allowed velocities (v x , v y ) and hence the al- 
lowed displacements, are not orthogonal to the instanta- 
neous position (x,y — ut) of the bob relative to the instan- 
taneous position of the support. It is easy to verify from 
the above equation that the allowed velocity (v x , v y ) is 
equal to the sum of a velocity vector (v x , v y — u) perpen- 
dicular to the relative position of the bob with respect 
to the point of suspension (x,y — ut), and the velocity of 
the support (0, u). 

v = v t + u 

The allowed displacements are vectors collinear to al- 
lowed velocities. A virtual displacement being the differ- 
ence of two allowed displacements, is a vector collinear to 
the difference of allowed velocities. Hence it is tangential 
to the instantaneous line of suspension. 

dr = vdt = v t dt + udt 

Sr = (v - v')dt = (v t - v' t )dt 

At any given instant string provides a tension along its 
length, with no shear (ideal string). Hence the constraint 
force, tension, still does zero work on virtual displace- 
ment. 

T ■ dr ^ 0, T • 5r = 

If one considers the support moving in a horizontal (or in 
any arbitrary direction), one can show that the allowed 
displacement is not normal to the instantaneous line of 
suspension. But the virtual displacement as defined in 
this article always remains perpendicular to the instan- 
taneous line of support. 

C. Motion along a fixed inclined plane 

The constraint is more conveniently expressed in the 
polar coordinate. The constraint equation is, 

f(r,0) = 0-0 o = 

where 0q is a constant. Hence the constraint equation for 
allowed velocities, Eq.(J2J, gives, 

Thus the allowed velocities are along the constant 
plane. Allowed velocity, allowed and virtual displace- 
ments are, 

v = fr, dr = rrdt, Sr = (f — r')rdt 



dr'= v'dt 
8r=(v - v') dt J/ 



dr =v dt^ 

/ i ° 



FIG. 3: Allowed and virtual displacements for a particle on a 
stationary slope 

If the inclined slope is frictionless (ideal), the constraint 
force provided by the surface is the normal reaction; 
which is perpendicular to the plane. Hence the work done 
by this force on allowed as well as virtual displacement 
is zero. 

N • dr = 0, N • Sr = 

D. Motion along a moving inclined plane 

For an inclined plane moving along the horizontal side, 
the constraint is given by, 

( x + ut ) ., a s n 
cot (6*0 ) = 

y 

f(x,y) = (x + ut) - cot(6 )y = 

whence the constraint for allowed velocities Eq.|J2J be- 
come, 

(x + u) — cot(0 o )y = 

Hence the allowed velocity (x,y) is the sum of two vec- 
tors, one along the plane (i + u, y), and the other equal 
to the velocity of the plane itself (— u, 0). 

V = Vj + u 

Allowed displacements are vectors along the allowed 
velocities, however the virtual displacement is still a vec- 
tor along the instantaneous position of the plane. 

dr = (v t + u)dt, dr' = (v' t + u)dt 
Sr = (v - v')dt = (v t - v' t )dt 

For the moving frictionless (ideal) slope, the constraint 
force provided by the surface is perpendicular to the 
plane. Hence the work done by the constraint force on 
virtual displacement is remains zero. 

N • dr y£ 0, N ■ Sr = 
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dr = v'dt 



dr=v dt ' 




v')dt 



u 

(velocity of 
the wedge) . 



FIG. 4: Allowed and virtual displacements for a particle on a 
moving slope 



IV. LAGRANGE'S METHOD OF 
UNDETERMINED MULTIPLIERS 

A constrained system of particles follow the equation 
of motion given by, 

m fc a fe = F fe + R fc , k = l,2,...,N 

where nik is the mass of the kth particle, a& is its accel- 
eration. F/j and R& are the total external force and force 
of constraint on the particle. If the constraints are ideal, 
we can write 



N 



J2 R* -Sr k =0 



(13) 



fc=i 



whence we obtain, 

N 



^(m fc a fe - F fe ) • Sr k = 



(14) 



k=l 



If the components of 5r k were independent, we could re- 
cover Newton's Law for unconstrained system from this 
equation. However for a constrained system 5rk are de- 
pendent through the constraint equations, 



/i(ri,r 2 , ...,r N ,t)=0, 



or. 



k=l 



dr k 



i = 1,2, 



i=l,2,. 



(15) 



(16) 



We multiply the above equations, Ea. (|16|) . successively 
by s scalar multipliers (Ai, A 2 , . . . A s ), called the La- 
grange's multipliers, and subtract them from the zero 
virtual work equation, Ea. H13|) . 



N 



E K-5> 



fc=i 



Sr k = 



(17) 



Explicitly in terms of components, 



N 

E 

k=l 



Ea. 



dx k 



Sx k + [Y} k Sy k + [Z] k Sz k j = 

( 18 ) 

where [Y) k and [Z] k denote the coefficients of Sy k and 
5z k respectively. 

The constraint equations Eq. fTB^I allows us to write the 
(s) dependent virtual displacements in terms of the re- 
maining n — 3N — s independent ones. We choose (s) 
multipliers (Ai, A 2 , . . . , A s ) such that the coefficients of 
(s) dependent components of virtual displacement van- 
ish. The remaining virtual displacements being indepen- 
dent, their coefficients must vanish as well. Thus it is pos- 
sible to choose (Ai, A 2 , . . . , A s ) such that all coefficients 
( [X] & , [K] & , [Z] & ) of virtual displacements (Sxk,Syk,Szk) in 
Ea. (|18|) vanish. Hence we have the forces of constraint 
in terms of the Lagrange's multipliers. 



E^#' 

1=1 



fc=l,2, 



,N 



(19) 



Thus the problem of mechanics reduces to finding solu- 
tion of equations of motion, 



m k a k =F k + ) Ai— , 
or k 

2—1 



~FT~~i fe = l,2, 



,N (20) 



with the constraints, 
/i(ri,r 2 , ...,r N ,t) 



0, 



1,2, 



, s 



(21) 



Thus we have to solve 3 AT + s scalar equations in 3iV + s 
unknown scalar quantities (x k ,y kl z kl Xi). After solving 
this system we can obtain the forces of constraint R& 
from Ea.fllty. 



V. GENERALIZED COORDINATES AND 
LAGRANGE'S EQUATIONS OF MOTION 

For the sake of completeness we discuss very briefly 
Lagrange's equations in generalized coordinates (for de- 
tail seo 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ). Consider a system of N 
particles under s holonomic, rheonomous constraints of 
the form given by Eq.JTJ. We can in principle express 
s of these coordinates in terms of the remaining 3N — s 
independent ones. Or we may express all the 3N scalar 
components of position in terms of n = 3N — s indepen- 
dent parameters q±, g 2 , . . . , q n and time (t). 

r k = r k ( qi ,q 2 ,...,q n ,t), k = l,2,...,N (22) 

The allowed and virtual displacements are given by, 

dr k 



dvi 



It a 

x dq~ Sqi ' 
3 = 1 J 

n „ 

EOY k s 
dq- 6 *' 

3 = 1 H] 



dt 



-dt, 



k = l,2,. 



, N (23) 
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From the Ea. lpHI) we obtain, 



VI. CONCLUSION 



(24) 

Introducing the expression of kinetic energy, 



1 N 



m k r 2 k 



fc=i 



and that of the generalized force, 



N 



dr k 



j = l,2,...,n (25) 



After some simple algebra one finds, 



^fddT dT , 



(26) 



since the qj are independent coordinates, coefficient of 
each Sqj must be zero separately. 



d dT _ dT_ 

dt dqj dqj 



j = l,2,...,n (27) 



In problems where forces are derivable from a scalar 
potential V(ri,r 2 , . . . ,rjv), 



In this article we make an attempt to present a quan- 
titative definition of the virtual displacement. We show 
that for certain simple cases the virtual displacement 
does zero work on forces of constraint. We also demon- 
strate that this zero work principle allows us to have a 
solvable class of problems. Hence we define this special 
class of constraint, the ideal constraint. We demonstrate 
in brief how one can solve a general mechanical problem 
by: i) Lagrange's method of undetermined multiplier and 
ii) Lagrange's equations in generalized coordinates. 

In Lagrange's method of undetermined multipliers we 
have to solve a larger number (3N + s) of equations, than 
in the case of Lagrange's equations (3 AT — s) in general- 
ized coordinates. However we can immediately derive the 
forces of (ideal) constraints in the former case. 

It is interesting to note that both the abovementioncd 
methods require the zero virtual work by constraint 
forces as a crucial starting point. In the case of La- 
grange's method of undetermined multipliers we start 
with the ideal constraint condition Ea. (|13|l . From there 
we write down Ea.lfH)l. Ea. fT7|l . Ea. l(T%|) and express 
the constraint forces in terms of Lagrange's multipliers, 
Ea. (|19|) . For Lagrange's equations in generalized coor- 
dinates we start with the ideal constraint, Ea. l(13|l . We 
work our way through Eq. (fH|l . Ea. (|2^|l . Eq. O^) and fi- 
nally obtain Lagrange's equations in generalized coordi- 
nates, Eq.(|23 and Eq.(|STJ- 



F fc = -V k V(T U T 2 ,...,T N ), k = l,2,...,N (28) 
we can write the generalized force as, 
'dr k \ dV 



Qj 



-v h v- 



dqj 



j = 1,2,..., n 
(29) 

Where V is the potential V expressed as a function 
of (gi, q2, ■ ■ ■ , q n ). In addition if the potential V does 
not depend on the generalized velocities, we obtain from 
Eq.©, 



d d(T - V) d(T - V) 



dt dqj 



dqj 



-0, j = l,2,...,n (30) 



At this stage one introduces the Lagrangian function L = 
T — V and in terms of the Lagrangian, the equations of 
motion Ea. 1)30(1 take up the form 



d dL dL 
dt dqj dqj 



0, j = l,2,...,n (31) 
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